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ABSTRACT 

We present limits on the cosmological constant from the statistics of gravitational 
lenses using newly completed quasar surveys, new lens data, and a range of lens models. 
The formal limit is Aq < 0.66 at 95% confidence in fiat cosmologies, including the 
statistical uncertainties in the number of lenses, galaxies, quasars, and the parameters 
relating galaxy luminosities to dynamical variables. The limit holds for either softened 
isothermal or de Vaucouleurs models of the lens galaxies, suggesting that the radial 
mass distribution in the lens galaxy is not a significant systematic uncertainty. The 
cosmological limits are unaffected by adding a small core radius to the isothermal lens 
models or by the luminosity of the lens galaxy. Inconsistent models of the effects of a 
core radius make significant errors in the cosmological constraints. Extinction in E/SO 
galaxies can significantly reduce the limits on the cosmological constant, but changing 
the expected number of lenses by a factor of two requires a minimum of 100 times 
more dust than is observed locally in E/SO galaxies. Rapid evolution of the lenses 
is the most promising means of evading these limits. However, a viable model must 
not only reduce the expected number of lenses but do so without significantly altering 
the average image separations, magnitudes, redshifts, and the relative properties of 
optically and radio selected lenses. 

Subject headings: cosmology: observations - galaxies: distances and redshifts - grav- 
itational lensing - galaxies: structure - galaxies: kinematics and dynamics 

1 Introduction 

Ostriker & Steinhardt (1995) invoke a "cosmic concordance" of models and data for flat 
cosmologies with a cosmological constant of Aq — 0.65 it 0.1. This value is already in dis- 
agreement with the most recent analysis of the statistics of gravitational lenses, which found 
that Aq ^ 0.5 at 90% confidence (Kochanek 1995). The statistical models were, however, 
limited in scope and included neither all the possible uncertainties nor recently discovered 
gravitational lenses or completed surveys, although the study treated some systematic uncer- 
tainties such as the mass distribution in E/SO galaxies. The current concordance is already 
a retreat from earlier invocations of Aq ~ 0.8 (e.g. Efstathiou, Sutherland, & Maddox 1990) 
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that were ruled out in the first full analyses of the observed statistics of gravitational lenses 
(Kochanek 1993, Maoz & Rix 1993). The dependence of the number of lenses on the cosmo- 
logical constant is steeper at Aq ~ 0.8 than Aq ~ 0.6, so it is possible that a careful treatment 
of the uncertainties will broaden the error bars sufficiently to permit the concordance. 

Turner (1990) first pointed out that cosmological models with a large cosmological 
constant produce large numbers of gravitational lenses. A fiat cosmology with Aq = 1 has 
roughly ten times as many gravitational lenses as a fiat cosmology with Aq = 0. A fiat model 
with Aq = 0.5 and an open model with Aq = and Oq = both predict about twice as many 
lenses as a fiat model with Aq = and Oq = 1- We only treat fiat models in this paper, 
but results for open cosmologies with Aq = can be approximated by mapping the range 
0.5 < Oo < 1 with Oo + Aq = 1 into the range < Oq < 1 with Aq = 0. The full dependence 
of lens statistics on both Aq and Oq is given in Carroll, Press, & Turner (1992) or Kochanek 
(1993). 

The first limits on the cosmological constant were derived from surveys for lensed quasars 
by Crampton et al. (1992), Yee et al. (1993), the Snapshot Survey (Bahcall et al. 1992, 
Maoz et al. 1992, 1993ab), and the ESO/Liege Survey (Surdej et al. 1993). These surveys 
examined 648 quasars and found five gravitational lenses. At the same time the MG (Burke 
et al. 1992) and JVAS (Patnaik et al. 1992) radio surveys found eight lenses. There are two 
in depth analyses of these data sets. Kochanek (1993) analyzed all the optical surveys and 
the radio lenses, but focused on a single lens model, and Maoz & Rix (1993) analyzed only 
the Snapshot Survey data while exploring a wide range of mass models and observational 
scaling laws. 

Kochanek (1993) used only the singular isothermal sphere (SIS) model for the lenses 
and focused on cosmological limits and dynamical normalizations. The upper limit on the 
cosmological constant was Aq < 0.8 at 90% confidence using systematic assumptions biased 
in favor of a high cosmological constant. The characteristic dark matter velocity dispersion 
of an galaxy was (Tdm* — 245 it 30 km s~^ at 90% confidence. This value was lower than 
predicted by Turner, Ostriker & Gott (1984), who estimated that the velocity dispersion of 
the dark matter should be (3/2)^/^ larger than the central velocity dispersion of the stars 
(cdm* — 275 km s~^). Better dynamical models by Kochanek (1993) and Kochanek (1994) 
(see also Breimer & Sanders (1993) and Franx (1993)) demonstrated that the assumptions 
leading to the (3/2)^^^ correction were incorrect; dynamically consistent models have (Tdm* — 
225 lb 10 km s~^, consistent with the estimates from fitting the observed lens separations. 

Maoz & Rix (1993) emphasized the effects of the mass distribution in the E/SO galax- 
ies rather than the cosmological limits, although they found a limit of Aq ^ 0.7 at 95% 
confidence. When Maoz & Rix (1993) tried de Vaucouleurs (1948) models for the lens 
galaxies using the dynamically estimated mass-to-light ratio of [M/ L)b — (10 it 2)h where 
Hq = 100/i km s~^ Mpc~^ (e-g- van der Marel 1991), they found that the models were in- 
capable of producing the observed image separations. Kochanek (1995) later demonstrated 
that [M I L)b — (22 ±3)h was required for de Vaucouleurs mass models to fit the lens sep- 
arations. Maoz & Rix (1993) then added a softened isothermal halo, and models with a 
sufficiently large halo could produce the observed image separations. The halo models were 
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not designed to be consistent with stellar dynamics (Breimer & Sanders 1993), but the total 
mass distribution of the successful models closely resembles the isothermal spheres used by 
Kochanek (1993). 

The dynamical conclusions of the statistical models {(Tdm* — Cc and [M/ L)b — 22h) are 
supported by other observations. Detailed models of the radio ring lens MG 1654+1346 by 
Kochanek (1995) also found that a nearly singular isothermal model fit the data extremely 
well and required a velocity dispersion comparable to the estimated dispersion of the stars, 
while a de Vaucouleurs model fit the data less well and required a mass-to-light ratio of 
[M/L)b — (20.4 lb 2.8)h. X-ray studies of ellipticals (e.g. Fabbiano 1989) and studies of 
disks and rings (e.g. Kent 1990; de Zeeuw 1991) also suggest that galaxies are embedded 
in dark matter halos, and stellar dynamical models are consistent with both constant mass- 
to-light ratios and dark matter (e.g. de Zeeuw & Franx 1991; van der Marel 1991; Saglia, 
Berlin, & StiaveUi 1992). 

Since these analyses, the NOT survey (Jaunsen et al. 1995) and the FKS survey 
(Kochanek, Falco, & Schild 1995) have increased the survey sample to 864 quasars with 
z > 1 and added one more probable lens (LBQS 1009-0252, Hewett et al. 1994, Surdej et 
al. 1994) to the statistical sample. At the same time, the MG (see Patnaik 1994), JVAS 
(see Patnaik 1994), and CLASS (Jackson et al. (1995), Meyers et al. (1995)) surveys have 
reached a total of eleven radio lenses. While we cannot accurately compute the absolute 
probability a radio source is lensed, we can reliably estimate the relative probability of the 
image separations. The observed properties of the radio lenses are powerful constraints on 
statistical models of the quasar surveys (Kochanek 1993). 

Several other systematic errors are treated in the literature. Galaxy evolution was 
treated by Mao (1991), Mao k Kochanek (1994) and Rix et al. (1994), who found that 
the conclusions of the no-evolution statistical models are correct unless galaxy evolution is 
still very dramatic at 2; ^ 1. Appealing to evolution to save models with a high cosmolog- 
ical constant is awkward, because one effect of a cosmological constant is to make galaxies 
form and evolve earlier than in a standard Oq = 1 CDM cosmology (see Carroll, Press, 
& Turner 1992). Moreover, most simple merger models preserved the expected number of 
lenses. Kochanek (1991) discussed the selection effects associated with quasar surveys, the 
brightness of the lens galaxy, and reddening. For lens surveys of bright quasars, we expect 
only reddening to be a serious concern. Some of these effects were reconsidered by Tomita 
(1995) and Fukugita & Peebles (1995). 

Kochanek (1992) noted that the relative redshifts of lens galaxies and their sources 
are an independent method of determining the cosmological model from the statistical ap- 
proaches used to analyze the quasar surveys. The test is weaker than a full analysis, but 
even it suggests that Aq ^ 0.9. Recently King (1994) (also Helbig & Kayser (1995)) strongly 
criticized this analysis for ignoring the measurability of the lens redshift and claim that the 
observable lens redshifts are insensitive to the cosmological model. Since the lens redshift 
test avoids some systematic errors in the standard statistical analyses (particularly problems 
associated with completeness, the number density of galaxies, and the average extinction) it 
is worth reevaluating the test in light of the critique. 
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The primary goal of this paper is to set firmer hmits on the cosmological constant in fiat 
cosmological models. In §2 we discuss the statistical model and its observational inputs. In 
§3 we discuss the effects of the uncertainties in the observational inputs on the cosmological 
limits using a softened isothermal sphere to model the lens galaxies. In §4 we reconsider the 
de Vaucouleurs lens models to see whether the cosmological limits are affected by changes 
in the radial structure of the lens galaxies. In §5 we continue our analysis of the softened 
isothermal models by examining the effects of finite core radii on both dynamical models 
of galaxies and lens statistics. In §6 we consider the effects of spiral galaxies, and in §7 we 
consider the effects of the lens galaxy luminosity and extinction. In §8 we reevaluate the 
lens redshift test, and in §9 we summarize the results. 



2 Methodology, Data, & Models 



In this section we outline the calculations and elaborate on improvements and changes 
from earlier discussions by Kochanek (1991, 1992, 1993) in the selection effects models, 
calculation techniques, statistics, number counts, and models for galaxies. We start by 
reviewing the numerical and statistical approaches. Next we discuss lens models and their 
parameters. Finally we discuss the number counts of galaxies and quasars. 



2.1 Calculation of Likelihoods 



If the two-dimensional lens potential is a function of radius 
is simply 

X d(j) 



u 



^ ^1 

r or 



then the lens equation 
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where u is the angular coordinate on the source plane and x is the angular coordinate on the 
image plane. The two-dimensional potential cf){r) satisfies the Poisson equation in angular 
coordinates V^^ = 2E/Ec where Ec = Dq^Dqi^I 4:TTGD^g is the critical surface density per 
square angle. The angular diameter distances -D05, Dqli ^^^^ ^LS relate proper distances 
to angular sizes between the observer and the source, the observer and the lens, and the 
lens and the source respectively (see Schneider et al. 1992). -Dql, Dls and Dqs denote the 
proper motion distances between the observer, lens, and source planes. 

If the surface density declines monotonically, the lens has two critical radii on which the 
magnification diverges (see Schneider, Ehlers, & Falco 1992, Blandford & Narayan 1986). 
The inverse magnification is 



1 



r dr 



1 



where the two factors are the eigenvalues of the inverse magnification tensor. The outer, or 
tangential, critical radius r_|_ is the solution of r — = 0, and the inner, or radial, critical 
radius r_ is the solution of 1 — (p^rr = 0. The caustics associated with the two critical lines 
are = for the tangential critical line (a pseudo-caustic technically) and m_ for the radial 
critical line. When two images are merging on the radial critical line, the third image is 
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located on the same side of the lens as the source at radius Vout- Sources inside the radial 
caustic are triply imaged, and sources outside the radial caustic are singly imaged. If we 
consider sources lying along an axis with u > 0, then the outermost image 1 lies on the 
same side of the lens as the source between the tangential critical line and the radius of 
the last multiple image (r_|_ < Xi < Vout)- The other two images are on the opposite side 
of the lens from the source, with image 2 lying between the radial and tangential critical 
lines (— r_|_ < ^2 < and image 3 lying between the origin and the radial critical 

line (— r_ < X3 < 0). When the source lies on the tangential caustic, images 1 and 2 are 
infinitely magnified creating an Einstein ring, and when the source lies on the radial caustic 
image 2 and 3 are infinitely magnified. An extensive discussion of the asymptotic statistical 
properties of these lenses is given in Blandford & Kochanek (1987). 

The comoving volume element at proper motion distance Dql is 

(1 + n,Dijrj,y/^ 

where the "curvature density" is 0/^- = 1 — Oq — ^0, the matter density is Oq, and the 
cosmological constant is Aq (Carroll, Press, & Turner 1992). If the total magnification of the 
three images for a source at position u is Af, then the probability, p[m^Zs)^ that a source 
with magnitude m and redshift z is lensed, assuming a survey selection function S{u) and a 
configuration selection function C(m), is 



/■°° , dn , /■"- udu ^dN , , , 

p(m, z,s) = 2Tr dL— dV uCu — m + 2.5 log M u 

Jo dL Jo Jo 47r dm 



dN 
dm 



(4) 



The comoving density of lenses labeled by luminosity L is dn/dL^ and the image positions, 
magnifications, and caustics must be computed for each luminosity and redshift. The location 
of the outer caustic (u-) and the magnification are functions of the lens model, the redshifts, 
and the luminosity. 

The selection function S{u) is equal to one when the current pattern of lensed images is 
detectable given their magnitudes, separations, and the survey selection function, and zero 
when it is not detectable. We consider the detectability of images in pairs (e.g. Webster et al. 
1988, Kochanek 1991), and an unresolved close pair (typically the 2-3 images) is merged into 
an average image before comparing it to the remaining image. The configuration selection 
function C{u) selects lenses with particular characteristics such as the image separations 
or the visibility of the central image. To estimate the probability of a particular image 
separation we calculated the probability of finding a lens within ±10% of the observed 
separation. 

The calculation is numerically simpler if we integrate over the image plane positions of 
image 2 rather than over the source positions. This helps stablize the integration because the 
Jacobian of the transformation from the image plane to the source plane introduces a factor 

into the integral and suppresses the near divergences of the integrand near the caustics. 
We also need only solve the lens equation (1) for the positions of the 1 and 3 images. The 
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lens probability with the change of integration variables is 

[rlN 1"^ 

(5) 



1 fin r^s rr+ rrlr rlN 



dN ^ 
dm 



where r is the coordinate of image 2, u is the source coordinate of the image, M2 is the 
magnification of image 2, and M is the total magnification of all three images. We selected 
the gridding in the three dimensions to resolve both the radial and tangential caustics and to 
keep the fractional errors smaller than 1%. For general lens models with a selection function, 
no part of this three dimensional integral simplifies, and it must be evaluated separately for 
each survey quasar. Most approximations found in the literature to simplify the calculation 
introduce quantitative errors comparable to the differences between Aq = and Aq = f/2 
cosmological models. 

2.2 Mass Distributions and Dynamical Normalization 

We use two mass models for the analysis: the de Vaucouleurs (f948) model, and a 
softened isothermal sphere model. The de Vaucouleurs model roughly matches the observed 
surface brightness of galaxies and it is the prototypical model for mass distributions with 
a constant mass-to-light ratio. The softened isothermal sphere is typical of dark matter 
dominated models and the inferred mass distribution of spiral galaxies. The hybrid models 
consisting of a de Vaucouleurs model combined with a softened isothermal sphere halo model 
used by Maoz & Rix (f993) and Breimer & Sanders (f993) so closely resemble a single 
softened isothermal sphere model that we do not treat them separately. 

The density profile of the softened isothermal sphere is 

a? 



'DM 



27rG'(r2 + s2 



(6) 



and it defiects light rays by 



d(p ^ ^ {r^ + s^f'"^ - s 
dr r 

where (Tdm is the one-dimensional velocity dispersion of the dark matter, s is the core radius, 
and b = ^tt{(Tdm I D^s I Dqs (Hinshaw & Krauss f987). This model was used by Fukugita 
& Turner (f99f), Fukugita et al. (f992), Krauss & White (f992), Maoz & Rix (f993), and 
Kassiola & Kovner (f993).^ The strength of the lens is determined by the ratio /3 = and 
the lens produces multiple images if /3 > f /2. The cross section of the model can be computed 
analytically (Hinshaw & Krauss f987), and under the assumption of a constant comoving 
core radius the cross section can be integrated to compute the optical depth (Krauss & White 
f992). 



The other softened isothermal model found in the literature is the two-dimensional lensing potential (f) = 
b{P + s^y/^ (Blandford & Kochanek 1987, Kochanek 1991, Wallington & Narayan 1993). Its lensing 
properties are virtually indistinguishable from the model of equations (6) and (7) if the Blandford and 
Kochanek (1987) core radius is twice the Hinshaw & Krauss (1987) core radius. 
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We assume the velocity dispersion (Tdm is related to the luminosity Z by a "Faber- 
Jackson" law of the form (Jdm = '^dm*{L/ L^Y^"' ^ in analogy to the Faber- Jackson (1976) 
law for the central velocity dispersions of F/SO galaxies. For a singular isothermal mass 
model, Kochanek (1994) found that (Tdm* = 225 it 10 km s~^ and 7 = 4.2 it 0.2, compared to 
(Tc* = 228ibl4 km s~^ and 7 = 3.4ib0.3 for the central velocity dispersions. We adopt broader 
uncertainties of (Tdm* = 225 it 22.5 km s~^ and 7 = 4.0 it 0.5 for the remainder of the paper 
to compensate for systematic errors in the dynamical models. The fact that (Tdm*I (^c* — 1 
demonstrates that the simple dynamical model with (Tdm = (3/2)^/^(Tc introduced by Turner, 
Ostriker, & Gott (1984) and used by Turner (1990), Fukugita & Turner (1991), Fukugita et 
al. (1992), Krauss & White (1992), Wallington & Narayan (1993), Maoz & Rix (1992), and 
Kassiola & Kovner (1993) is incorrect. Using the Turner, Ostriker, & Gott (1984) dynamical 
model overestimates the number of lenses by 125% and overestimates their separations by 
50%. 

It is increasingly clear that F/SO galaxies are effectively singular. HST observations show 
cores with surface brightness distributions approaching with ^ j/ ^ 1 (see Tremaine et 
al. 1994). This range matches that between the Hernquist (1990) and Jaffe (1983) models, 
where the Jaffe (1983) model is the same as a singular isothermal sphere in its inner regions.^ 
The almost uniform absence of central images in the observed lenses (Wallington & Narayan 
1993, Kassiola & Kovner 1993) and models of individual lenses (e.g. Kochanek 1995) imply 
core radii smaller than 100/i~^ pc for galaxies."^ This suggests that earlier concerns about 
the effects of finite core sizes on lens statistics (e.g. Blandford & Kochanek 1987, Kochanek 
k Blandford 1987, Hinshaw k Krauss 1987, Kochanek 1991, Krauss k White 1992) were 
overstated. Nonetheless, we consider models with core radii, where is the core radius for 
an galaxy, and s/s* = (cr/cr*)^ with x = 4.8. For a "Faber- Jackson" exponent 7 = 4 
this corresponds to s/s* oc [L/ L^Y''^ and matches the observed scaling of the effective radius 
with luminosity (Fukugita k Turner 1991). 

The second model we consider is the de Vaucouleurs (1948) model, as an example of 
a constant mass-to-light ratio model. The lens defiection produced by the de Vaucouleurs 
model is 

— = g[R Re) where g[x) = , , and k = 7.67. 8 

or c^Ke Uos X Jq ax X exp{ — kx'^/^) 

The characteristic defiection scale is 

4Gm Dls 



3'.'9 



OS 



m 



h ^kpc 

Rp 



OS 



for total mass m and effective radius Re. The defiection function peaks at Rmax = 0.164i?, 
with max(^,r ) = 0.746^,. The effective radius of F/SO galaxies varies with luminosity as Re = 



^The core radius s of the dark matter need not match the core radius of the luminous matter, but there are 
dynamical problems if the two differ by any large factor. 

^Steeply dropping radial profiles like the Hubble profile will have larger, finite core radii, but similar high 
central surface densities, (e.g. Nair, Narasimha, & Rao 1993). 
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Re*{L/L^Y with i?e* = (4 ± l)h-^ kpc and a = 1.2 ± 0.2 where Hq = lOOh km s'^ Mpc'^ 
(e.g. Kormendy & Djorgovski 1989, Rix 1991). The average blue mass-to-light ratio varies as 
{M/L)b = {M/L)b*{L/L^)\ where {M/L)b* = (10 ± 2)h and b = 0.25 ± 0.10 in dynamical 
models assuming constant mass-to-light ratios (van der Marel 1991, Kormendy & Djorgovski 
1989). With M^{B^) = -19.9;o:4 (Efstathiou, EUis, & Peterson 1988, Loveday et al. 1992) 
the deflection scale becomes K = (1''3 ± 0'.'4:){L/ L^ f-°^^°-^^DLs/ Dqs- The magnitude of the 
deflection scale compared to the observed lens separations is the fundamental reason that 
Maoz & Rix (1993) and Kochanek (1995) found that the de Vaucouleurs models could not 
flt the observed image separations. 

Maoz & Rix (1993) and Breimer & Sanders (1993) also examined hybrid models by 
adding an isothermal halo with a large core radius to the de Vaucouleurs model. When 
these models are examined in detail, the hybrid model so closely mimics a softened but 
nearly singular isothermal sphere that there is little point in making the artiflcial division 
into a luminous mass distribution and a dark matter halo with the resulting multiplicity of 
parameters. In the hybrid models the halo determines the asymptotic deflection and the 
inner de Vaucouleurs model determines the effective core radius. For example, the deflection 
produced by a de Vaucouleurs model inside the peak at Rmax = 0.164i?e matches that of 
an isothermal sphere with asymptotic deflection equal to the peak deflection and a core 
radius of s = 0.0072i?e = 30/i~^(iv/iv*)^'^ parsecs. With the addition of a halo the dynamical 
normalizations must be recomputed (this was done by Breimer & Sanders (1993) and Franx 
(1993) but not by Maoz & Rix (1993)) with the unsurprising result that the normalizations 
match those derived for isothermal spheres by Kochanek (1993, 1994). 

We want to use lens models that are dynamically consistent, or at least to know when 
the models capable of fltting the lens data are inconsistent with the requirements of stellar 
dynamics. Where the dynamical normalization is unknown, we flt constant isotropy dynami- 
cal models to the sample of 37 F/SO galaxies from van der Marel (1991). If p is the luminosity 
density of the galaxies found by deprojecting the observed surface density, /3 = 1 — (tH (t^ is 
the constant isotropy, and M[r) is the mass, then the radial velocity dispersion is 



where z is the coordinate along the line of sight, R is the projected radius, and E is the 
observed surface brightness of the galaxy (Binney & Tremaine 1987). For the purposes of 
this paper we flt isotropic models (/3 = 0) but broaden the errors on the dynamical variables 
by approximately a factor of two. 



■*Note the typographical errors in equations (2) and (3) of Kochanek (1994). Equation (2) should read va^ 
and equation (3) should read • 




(10) 



and the line of sight velocity dispersion is 
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2.3 Statistical Analysis 



Given a lens model, the number counts of galaxies, the properties of the galaxies, and 
the lens survey selection function, we compute the probability pi that quasar i is lensed or 
the probability pi[A6i) that quasar i is lensed and has image separation A^^ using equation 
(5). If the sample contains Nu unlensed quasars, Nl lensed optical quasars, and Nr radio 
lenses, the likelihood function for the data d given a fixed set of model parameters ^ is 



We expanded the logarithm of the probability that a quasar is unlensed ln(l — pi) ~ —pi for 
Pi <^ 1, and we include only the relative probability that a radio lens has separation A^^ 
rather than the absolute probability that it was lensed by using the ratio pk{A6k) /pk- The 
various probabilities include all the model parameters and survey selection function models. 
We can add other aspects of the configuration of the lenses to the likelihood such as the 
lens galaxy redshift, the detectability of odd/central images, and the image morphology if 
desired. 

The model parameters also have uncertainties that are not included in any of the existing 
cosmological limits. For most parameters we have estimates of their mean values and disper- 
sions, so we can define the prior probability P{C]) for parameter ^j. The prior probabilities 
are assumed to be normal distributions for signed variables and log-normal distributions for 
positive definite variables. Using Bayes theorem, the probability of the parameters given the 
data is then 



The denominator of the expression merely normalizes the total probability to be unity. We 
use the Bayesian analysis because it provides a clear, formal approach for including the 
parameter uncertainties, for defining error bars on non-Gaussian distributions, and for pro- 
jecting out nuisance variables. The probability distribution for any subset of the parameters 
is found by integrating the probability distribution over the unwanted variables. 

We usually reduce the distribution to a two dimensional parameter space and present 
likelihood ratio contours in this subspace. Usually we examine the space of the cosmological 
model and the dynamical variable. The 68%, 90%, 95%, and 99% confidence limits on one 
parameter are the extrema of the likelihood ratio contours at 61%, 26%, 14% and 3.6% 
of the peak likelihood. The Bayesian limit on one variable is found by marginalizing the 
distribution to a single variable and finding the limits enclosing the relevant fraction of the 
probability. 






(13) 
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2.4 Number Counts of Galaxies 

We divide the number counts of galaxies into an E/SO class and spirals, with comoving 
densities fitted by a Schechter (1976) function, 

We use an average comoving density of galaxies of = (1.40 it 0.n)h^ Mpc""^ (Loveday et 
al. 1992) where Hq = lOOh km s~^ Mpc~^. The division of the luminosity function by galaxy 
type is taken from Marzke et al. (1994), renormalized to the total galaxy density found by 
Loveday et al. (1992). The Loveday et al. (1992) sample is too deep to type the fainter 
ellipticals (Marzke et al. (1994)), while the total density in the Marzke et al. (1994) sample 
is too high because of local structures. With these assumptions the comoving densities of 
E/SO and spiral galaxies are 

ne = (0.61 ± 0.21)/^^ X 10"^ Mpc-3 and = (0.79 ± 0.26)/i^ x 10"^ Mpc'^ (14) 

We adopt a common slope a = —1.00 it 0.15 for both galaxy types (Marzke et al. 1994), 
which also agrees with the total slope found by Loveday et al. (1992). These estimates 
supersede the (statistically and systematically) less accurate model of Eukugita & Turner 
(1991) based on the earlier number counts of Efstathiou, Ellis & Peterson (1988) and the 
morphological type ratios of Postman & Geller (1984). While the error bars in the newer 
surveys are more reliable, they are also broader. The quoted error estimates are all one 
standard deviation errors. The shortcoming of this error model is that the values of a and 
must be correlated, but neither Marzke et al. (1994) nor Loveday et al. (1992) discuss the 
parameter uncertainties in the full three-dimensional space of a, and L^. The absolute blue 
magnitude of an E/SO galaxy is M^[Bj) = —19.9^0^4 magnitudes, which corresponds to 
luminosity = lA^oi X 10^° Lb® (Efstathiou, EUis, & Pet erson 1988, Loveday et al. 1994). 
When we calculate the selection effects from the galaxy luminosity (§7.1) and the probability 
of measuring a lens galaxy redshift (§8), we use the A;-corrections model of Coleman, Wu, & 
Weedman (1980) to compute the apparent magnitude of the lens galaxy. The model does not 
include evolutionary corrections, and it is probably accurate only to within one magnitude. 

2.5 Number Counts of Quasars 

The number counts model for the quasars is critical to the statistical analysis because 
magnification bias (Gott & Gunn 1974, Turner 1980) has an enormous effect on lens prob- 
abilities in bright quasar surveys. No realistic calculation of lens probabilities can neglect 
its effects. In this section we compare the various number counts models used in lensing 
calculations. 

Turner (1990), Eukugita & Turner (1991), and Kochanek (1991, 1993) used a simple 
two power law model with a fixed apparent break magnitude developed by Turner (1990) 
based on the compilation of data by Hartwick & Schade (1990). This model has a break 
at niQ = 19.15 B mags, a bright end slope a = 0.86 and a faint end slope /3 = 0.28 where 
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log;^Q dN I dm = a[m — mo) + c for m < mo and log^^Q dN j dm = fi{m — mo) + c for m > mo- 
We are uninterested in the value of the constant c for lens calculations. The bright slope 
a in the Turner (1990) model is the slope for the apparent magnitude number counts of all 
bright, UV-excess selected quasars with z <^ 2.2 (Hartwick & Schade 1990). By averaging 
over the redshifts of the quasars, this model significantly underestimates the true slope of 
the luminosity function. 

Wallington & Narayan (1992) and Kassiola & Kovner (1993) use a quasar luminosity 
function modeled on that of Boyle et al. (1988). In this model the luminosity function is 

^ oc [lO-«("-"«) + iQ-Pi^-^o)] (15) 
dm I- J 

where the bright end slope is a = 1.12 it 0.06, and the faint end slope is /3 = 0.18 it 0.08. 
The break magnitude mo evolves with redshift based on the Boyle et al. (1988) model for 
redshifts less than three, and with a correction of —0.54(2; — 3) magnitudes above a redshift 
of three to fit the higher redshift surveys (Wallington & Narayan 1993).^ The Snapshot 
Survey and Maoz & Rix (1993) used a two power law model in absolute magnitude based 
on the luminosity function of Boyle et al. (1987). The bright end slope is a = 1.04, the faint 
end slope is /3 = 0.08, and the evolutionary model for mo is identical to that of Wallington 
& Narayan (1993). 

The complexity of the absolute magnitude models hides the fundamental simplicity 
of the evolution in apparent magnitude. Between the redshifts of 1 and 3, the apparent 
magnitude of the break is ~ 19 B magnitudes. Below a redshift of one and above three we 
allow the break magnitude to evolve as 

r mo + (^ - 1.0) z < 1 

mo{z) = <mo 1<2;<3 (16) 

[ mo - 0.7{z - 3) 3 < z 

This evolution model matches the Wallington & Narayan (1993), Boyle et al. (1988, 1990), 
and Maoz & Rix (1993) models to higher accuracy than required by the statistical un- 
certainties. We fit this evolution model and the broken power law functional form to the 
quasar luminosity function data in Hartwick & Schade (1990) for 2; > 1 to estimate the three 
parameters a, /3, and mo. 

Figure 1 shows the projected likelihood functions for the three parameters, with crosses 
at the Fukugita & Turner (1991), Maoz & Rix (1993), and Wallington & Narayan (1993) 
models. The best models with the simple apparent magnitude evolution model in equation 
(16) are slightly better than the more complicated evolution models used by Wallington & 
Narayan (1993) and Maoz & Rix (1993), although the improvements are not significant. 



®It appears that a conversion between Hq = 50 km s~^ Mpc~^ used by Boyle et al. (1988) and Hq = 
100 km s~^ Mpc~^ used by Wallington & Narayan (1993) was lost in their conversion to apparent magnitudes. 
This implies that figures as a function of apparent magnitude in Wallington & Narayan (1993) and Kassiola 
& Kovner (1993) (who used the Wallington & Narayan (1993) model) must be shifted by 1.5 magnitudes. 
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Fig. 1. -Likelihood contours for the quasar number counts model. In each two dimensional subspace 
the third variable was optimized. The three parameters are the bright slope a, the faint slope /3, 
and the break magnitude nio. Contours are shown at the 68%, 90%, 95%, and 99% confidence 
limits for one variable in the likelihood ratio. The crosses mark the best fit, the Fukugita & Turner 
(1991) model (FT), the Wallington & Narayan (1993)/Boyle et al. (1988) model (WN), and the 
Maoz & Rix (1993)/Boyle et al. (1987) model (MR). The FT model has too fiat a bright slope to 
be consistent with the data, the MR model has too fiat a faint slope to be consistent with the data, 
and the WN model is consistent with the data. 

The best model has a = 1.07 ± 0.07, /3 = 0.27 ± 0.07, and nio = 18.92 ± 0.16 B mags (90% 
confidence error bars) with = 63.6 for N^of = 62 degrees of freedom. A one standard 
deviation change in the from = N^of is Ax^ = 11 for 3 parameters. The bright end 
slope of a = 0.86 in the Fukugita & Turner (1991) model is in complete disagreement with 
the data - the likelihood ratio is log[L / L^ax) = —12.1. The best model with the Wallington 
& Narayan (1993) slopes has mo = 18.95 and log[L / L^ax) = —0.95. This is within the 
90% confidence interval, and our best fit is no better than the full Wallington & Narayan 
(1993) model using the Boyle (1988) functional form rather than the broken power law. The 
best model with the Maoz & Rix (1993)/Boyle et al. (1987) slopes has mo = 19.18 and 
\og{L I Lraax) = —2.5. It cau be rejected with 99% confidence because the faint slope of 
fi = 0.08 is inconsistent with the data and with later analyses (Boyle et al. 1988, 1990). We 
will use our best fit model (a = 1.07, fi = 0.27, and mo = 18.92 B mags) and the evolution 
model for mo in equation (16). We let the quasar density be constant beyond nic = 30 B 
magnitudes, and we use an average B — V color of 0.2 mag as suggested by Maoz et al. 
(1993ab). 
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Table 1: 


Summary of Lens Data 






Lens 




^/ 








A6' 


Comments 








mag 










1208+1011 


3 80 




1 7 9 




ZjVJ . t: 


0.45 




H1413+117 


2 55 




1 7 




1 9 5 

X tJ . 


1.2 




LBOS 1009-0252 


2 7^ 




1 7 9 




20 

ZjVJ . VJ 


1.6 




PQ 1115+080 


1 72 


29 


1 6 1 

X VJ . X 


X tJ . 


21.2 


2.2 




0142-100 


2.72 


49 

VJ . Tit./ 


1 6 8 

X VJ . 


1 9 

X tJ . VJ 


21.0 


2.2 




0957+561 


1.41 


0.36 


16.7 


18.5 




6.1 


Not Used 


B 0218+356 




68 




20 




0.35 




MG 0751 +271 fi 








1 Q 




9 

VJ . >J 






1 fi9 




1 ^ 




1 q n 


1 3 

X . ■ } 




CLASS 1600+434 


1.61 




20.0 




19.0 


1.4 




MG 1549+3047 




0.11 


23 


16 




1.7 




B 1938+666 






(23) 






1.8 




MG 1654+1346 


1.74 


0.25 




19.0 


21.0 


2.1 




MG 0414+0534 


2.64 




21.8 


r=2i.4 




2.1 




MG 1131+045 








(22) 




2.1 




CLASS 1608+656 


2.30 


(0.63) 


(20) 


(20) 


21.0 


2.1 




2016+112 


3.27 


(1.01), 


? i=21.8 


i=21.9 




3.8 


Not Used 



NOTES: The source and lens redshifts are and 2;;, the source and lens magnitudes are 
and m;, and the image separation or the diameter of the tangential critical hne is A^. The 
magnitude niio is our estimate of how bright the lens galaxy has to be for the lens redshift to 
be measurable for the systems used in §8. The first six lenses are in the optical surveys, and 
the last eleven lenses are in the radio surveys. An entry in parenthesis means low accuracy 
or substantial uncertainty. The magnitudes are R magnitudes unless otherwise specified. 
See Surdej & Soucail (1994) for more details and the original references. 

2.6 Observational Data & Selection Functions 

We use the optical lens surveys by Crampton et al. (1992), Yee et al. (1992), the 
Snapshot Survey (Bahcall et al. 1992, Maoz et al. 1992, 1993ab), the ESO/Liege Survey 
(Surdej et al. 1993), the NOT survey (Jaunsen et al. 1995), the HST GTO survey (see Ealco 
1994), and the EKS Survey (Kochanek et al. 1995) but include only the quasars with z > 1. 
We use the conservative evaluation of the ESO/Liege Survey candidates by Kochanek (1993), 
and the preliminary EKS selection function. The combined sample contains 862 quasars and 
six lenses. Table 1 summarizes the observational data. We do not use 0957+561 in the 
analysis because the lens consists of both a cluster and a galaxy. Kochanek (1995) and 
Wambsganss et al. (1995) examine the statistics of wide separation lenses like 0957+561. 
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In addition to the optical samples we include the image separations of the radio lenses 
found in the MIT/Greenbank (MG) Survey (see Burke et al 1992, Patnaik 1994), the Jodrell 
Bank Survey (see Patnaik et al. 1992, Patnaik 1994), and the CLASS survey (Jackson 
et al. 1995, Myers et al. 1995). Because the redshift and luminosity functions of these 
radio surveys are incompletely known, it is impossible to calculate the absolute probability 
of gravitational lensing in these samples reliably. We can, however, estimate the relative 
likelihoods of finding the observed image separations with much weaker uncertainties from 
the effects of the magnification bias on the likelihood of finding a particular separation. 
There are eleven of these radio lenses and Table 1 summarizes the observational data. We 
used all these lenses except 2016+112 in the analysis. The lens 2016+112 is probably a 
complicated hybrid lens involving two galaxies at different redshifts. Where the source 
redshift is unknown, it is set to = 2. The optical magnitudes are set to m = 19 for the 
magnification bias calculation. The bias near m = 19 is comparable to that expected for the 
radio samples. As discussed by Kochanek (1993) the source redshift has little effect on the 
image separations, particularly in fiat cosmological models. 

3 Isothermal Lenses and Parameter Uncertainties 

Independent of any systematic problems in the statistical model, we must consider how 
the uncertainties in the model parameters affect the cosmological conclusions. The known 
statistical uncertainties in galaxy number counts, quasar number counts, and dynamical 
relations have not been included in any of the full analyses of lens statistics. Kochanek (1993) 
included the uncertainties in the dynamical normalization for the lens galaxies and examined 
the effects of the other variables in a particular cosmological model, while Maoz & Rix (1993) 
examined "conspiracy" models in which parameters were biased to increase the likelihood of 
a large cosmological constant. In this section we derive the cosmological limits including the 
uncertainties of the input parameters. We use the only known lens model consistent with 
lens statistics, galaxy dynamics, and individual lens models, the nearly singular isothermal 
sphere, with a core radius of = 10/i~^ pc. This model also mimics the Breimer & Sanders 
(1993), Franx (1993), and Maoz & Rix (1993) models of an E/SO galaxy embedded in a dark 
matter halo. All later calculations will automatically include these uncertainties. 

3.1 Galaxy Number Counts and Dynamical Relations 

The variables relating galaxy number counts and the isothermal lens model are the 
number density of E/SO galaxies He = (6.1 ± 2.1)h^ X 10"^ Mpc~^ , the Schechter (1976) 
function slope a = —1.0 ± 0.15, the Eaber- Jackson (1976) exponent 7 = 4.0 ± 0.5, and 
the velocity dispersion of the dark matter for an galaxy, (Tdm* = 225 ± 22.5 km s~^ 
(Kochanek 1994). The prior probability distributions for the variables are assumed to be 
the log-normal distribution for the galaxy number density and the normal distribution for 
the other variables with the stated one-standard deviation uncertainties. If (-D05) 
sample average of the cube of the proper motion distance to the quasars, the expected 
number of lenses varies as ne{DQg)crjjj^^T[l + a + 4/7] and the mean image separation 
varies as crjjj^^T[l + a + 6/7]/L[l + a + 4/7] (see Kochanek 1993). Eigure 2 shows the two- 
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Fig. 2. -Likelihood contours showing the dependence of the cosmological limits on the parameters 
determining the galaxy number counts and dynamical relations. Contours are shown at the 68%, 
90%, 95%, and 99% confidence limits for one variable in the likelihood ratio. The dashed lines 
show the prior probability distribution for the variable (a, 7, We, or (Jdm*), the light solid lines 
show the limits using only the lensing data (although including the prior probability distributions 
for the other three variables), and the heavy solid lines show the joint limits including the lensing 
data and the prior probability distribution. The galaxy density is in units of h'^ Mpc""^. 
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dimensional likelihood contours for each of the four variables and the matter density Oq in 
a flat universe after marginalizing the original five dimensional likelihood function over the 
other three variables. 

With the inclusion of all the statistical uncertainties (artificially broadened uncertainties 
for 7 and (Tdm*) the two standard deviation upper limit on the cosmological constant is 
Ao < 0.66. The best fit models with Oq — 1 are consistent with the external estimates 
of the input parameters. In particular, the dynamical estimates of (Tdm* from Kochanek 
(1994) exactly match the best estimates from fitting the lens separations unless Aq ~ 1. 
The existence of a (3/2)^/^ correction for the dark matter velocity dispersion is supported 
neither by galaxy dynamics nor the observed separations of gravitational lenses. The best 
fit estimate is (Tdm* = 220 it 20 km s~^, and it is well constrained because the average image 
separation is a strong function of the velocity dispersion, oc crj^^^- 

The lensing data is least able to differentiate between changes in the galaxy number 
density Ue and the cosmological model. Because the image separation distribution for the 
singular isothermal sphere does not depend on cosmology, only the variation of the lens 
probability with redshift differentiates between the galaxy number density and cosmology. 
With only five optical lenses (the radio lenses do not matter here) there is too little data 
to measure the variation with redshift. This is a fundamental uncertainty, and any other 
systematic problem that can mimic changes in Ue is similarly difficult to constrain. For 
example, using dust to make some fraction of the E/SO galaxies completely opaque at all 
redshifts is the same as changing Ue- 

The Schechter (1976) function slope a controls the relative numbers of low and high 
mass galaxies. The lens statistics are dominated by the number of ^ galaxies as long as 
a ^ — 1, and in this regime the results will be insensitive to the exact value of the exponent. 
When a <^ —1 the number of lenses produced by the small galaxies rises rapidly, so there are 
strong restrictions on a ^ — 1. In reality there is a correlation between the value of Ue and 
the value of a. Unfortunately the studies of galaxy number densities (Marzke et al. 1994, 
Loveday et al. 1992, Efstathiou et al. 1988) give only the value of Ue and its uncertainties 
for the best fitting value of a, so we are unable to include it in the error model. 

The "Faber- Jackson" exponent is weakly constrained by the lensing data to lie near 
the standard values of 7 ~ 4 with strong constraints on 7 < 4. The expected number of 
lenses is minimized for 7 ~ 2.8 (for a = —1), and there the likelihood contours reach their 
peak values of the cosmological constant. The exponent 7 and (Tdm* are correlated, but 
we chose to use broad uncorrelated error bars for the two variables separately rather than a 
detailed correlation matrix that was dominated by systematic uncertainties in the dynamical 
modeling. 

3.2 Quasar Number Counts 

No study of the lens surveys has examined the sensitivity of the cosmological results to 
the quasar number counts model. Kochanek (1993) examined how well the lensing model 
could estimate the values of the individual parameters of the number counts model for a fixed 
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Fig. 3. -Likelihood contours showing the dependence of the cosmological limits on the quasar 
number counts model. Contours are shown at the 68%, 90%, 95%, and 99% confidence limits for 
one variable in the likelihood ratio. In the left panel the dashed contours are for the Fukugita & 
Turner (1991) number counts model and the solid contours are for the Maoz & Rix (1993)/Boyle 
et al. (1987) number counts model. In the right panel, the dashed contours are for our best fit 
number counts model, and the solid contours are an average of the best fit model plus eleven other 
models randomly selected from the probability distribution of the fits to the quasar number counts. 
The likelihood distributions include the uncertainties in (Jdm*, a, and 7. 



cosmological model and found that the lensing data was better fit by a steeper slope for the 
number counts of bright quasars than in the Fukugita & Turner (1991) model. Since there 
is no convenient way of storing intermediate results to allow a multi-dimensional parameter 
study as used for the distribution of lens galaxies, we restricted ourselves to a small Monte 
Carlo study. We selected 12 sets of number counts parameters from our fits to the Hartwick 
& Schade (1990) compilation, as well as the Fukugita & Turner (1991) number counts model 
and the Maoz & Rix (1993)/Boyle et al. (1987) number counts model, for a total of 14 
different models. The probability of fitting the separations of the radio lenses was kept fixed 
to avoid introducing any biases in the models. 

Figure 3 shows the likelihood contours in the space of (Tdm* and Oq including the prior 
probability distributions for ne, a, 7, and (Tdm*- The best fit model derived in §2.5, the Maoz 
& Rix (1993) model, and the Wallington & Narayan (1993) model all produce similar results, 
and averaging over the small Monte Carlo sampling of the distribution has little effect on the 
cosmological conclusions. The Fukugita & Turner (1991) model, however, produces different 
cosmological limits because of its fiatter bright quasar slope. The Fukugita & Turner (1991) 
number counts model give 68%, 90%, 95% and 99% confidence limits on the cosmological 
constant of Aq ^ 0.58, 0.75, 0.84, and 0.90 respectively while the other three models give 
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Fig. 4. -The left panel shows the 95% confidence contour for the the mass and effective radius 
i?e* of an galaxy for JIq = 1 (light solid) and Aq = 1 (light dashed) when a' = 1. The dashed 
(solid) heavy solid lines show the 95% confidence contours for JIq = 1 and a! = 0.8 (a' = 1.2). 
The likelihoods include the observed log-normal prior probability distribution for i?e* and they are 
drawn relative to the peak for each model. The cross shows the observational uncertainties for 
Re* and m* in constant mass-to-light ratio models. The right panel shows the likelihood contours 
for the mass of the galaxy and the cosmological model. The uncertainties in the galaxy number 
counts (a and We), the normalization of the effective radius i?e*, the mass scale m* and the scaling 
exponent of the mass-to-light ratio b are included in the error bars. The light solid lines show 
the 68%, 90% , 95%, and 99% confidence intervals for one parameter in the likelihood ratio when 
a' = 1, and the dashed (solid) heavy lines show the 95% confidence interval for a' = 0.8 (a' = 1.2). 



Ao ^ 0.33, 0.57, 0.66, and 0.78 respectively. The underestimate of the magnification bias 
by the Fukugita & Turner (1991) number counts model is the reason that Kochanek (1993) 
found weaker limits on the cosmological constant than Maoz & Rix (1993). 

Because the spread in the maximum likelihoods is only a factor of two, we cannot 
choose between the luminosity function models. The sample has one bright (PG 1115+080), 
one faint (1208+101), and three intermediate (0142-100, LBQS 1009-0252, H 1413+117) 
magnitude quasars. The likelihoods are similar because the gains (or losses) in the probability 
of PG 1115+080 are compensated by a losses (or gains) in the probability of 1208+101. Thus 
the probability for the parameters of the quasar luminosity function are determined by the 
prior probability distribution for fitting the observed quasar number counts. We use only 
the best fit model from §2.5 for the remainder of the paper. 
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4 De Vaucouleurs Models 



From both lens models (Kochanek 1995) and lens statistics (Maoz & Rix 1993, Kochanek 
1995) we know that the de Vaucouleurs mass model for lens galaxies with standard mass- 
to-light ratios is unable to fit lens data. One solution is to add a halo to the de Vaucouleurs 
model (Maoz & Rix 1993, see also Breimer & Sanders 1993), but these galaxy -|- halo models 
so closely resemble softened isothermal sphere models that we will not treat them separately. 
In this section we explore the properties of the de Vaucouleurs models to estimate the galaxy 
masses required to fit the lens data, and to see if the radically different shape of the de 
Vaucouleurs model has any effect on the cosmological limits. 

The model depends on seven parameters: a, ne, i?e, {Mj L)^, 6, and M^[Bt) (see 
§2.2). We combine the mass-to-light ratio (M/L)^ with the luminosity scale M^[Bt) and 
fit the mass of an galaxy m* as the fundamental parameter. We can do this because we 
do not have the covariance matrix for a, Ue and in the galaxy number counts models. 
The mass scales with the luminosity as m/m^ = (L/L^Y'^'^ where b is the exponent for 
the scaling of the mass-to-light ratio (M/L) = (M/ L)^{L/ L^Y . We tabulated the lensing 
probabilities as a function of mass and effective radius, which makes it easy to examine 
scaling laws of the effective radius with mass [Re oc with a' = a / [1 -\- b)) but hard to 
examine ones with luminosity [Re oc Z"). For our standard models we leave the mass m* as a 
free parameter to be determined from the lens data, and average over the uncertainties in the 
galaxy number counts model [a and ne), the effective radius of an galaxy [Re* = [4:±l)h~^ 
kpc), and the exponent of the mass-to-light ratio [b = 0.25 it 0.10) for a fixed value of 
a' = a/[l -\- b) = 0.96 it 0.25. We examine a' = 0.8, 1.0, and 1.2, which spans the range for 
a. The average over the value of i?e* can be viewed either as including the uncertainties in 
i?e* or as an average over the fundamental plane. 

The curious property of the standard scalings for the de Vaucouleurs model is that the 
peak defiection is proportional to 1}+^-"^ = ^o.osio.s compared to L^'^ in the isothermal 
models. The peak defiection is nearly the same for all luminosities, and either an L* galaxy 
can generate the largest observed separations or no galaxy can. The peak defiection is 
proportional to mj Re and by the virial theorem we expect that mj Re oc . Observed 
velocity dispersions have oc L^'^ so the standard de Vaucouleurs scalings do not match the 
Faber- Jackson (1976) law because either the slope of the mass-to-light ratio is too shallow, 
the slope of the effective radius scaling is too steep, or the measured velocity dispersion 
differs systematically from the virial dispersion. 

If we completely ignore the observed effective radii of galaxies and determine what de 
Vaucouleurs profile would (in abstract) best fit the data, it is a compact [Re ~ 400/i~^ 
pc), low mass (m* ~ 8 X l{)^^h~^ Mq) galaxy with Aq ^ 0.75 at 95% confidence. When we 
add the prior probability distribution for the observed effective radii, the best fit models for 
Oo = 1 are i?e* = 3.0/i~^ kpc and m* = 2.5 X 10^^ /^"^Af© for a' = 0.8, i?e* = 3.3/i~^ kpc and 
= 3.2 X W^h-^M^ for a' = 1.0, and i?e* = 4.0/^"^ kpc and = 4.0 x W^h-^M^ for 
a' = 1.2. As a' is lowered, the effective radius and mass decrease. The left panel of Figure 
4 shows the 95% confidence regions for the parameters. The models with Aq = 1 do not fit 
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the data nearly as well because the models are driven to low masses and small effective radii 
to maintain constant image separations while minimizing the number of lenses. The right 
panel of Figure 4 shows the cosmological limits for a! = 0.8, 1.0, and 1.2. In each case the 
best fit model is Oq = 1, and the 95% confidence upper limits on Aq are Aq ^ 0.68, Aq ^ 0.62, 
and Aq ^ 0.55. The results for a weighted average of the results for different a' are similar 
to those for a' = 1.0. 

In short, we confirm the results of Maoz & Rix (1993) using a larger data set and a 
more thorough statistical analysis. Constant mass-to-light ratio models of E/SO galaxies 
are inconsistent with gravitational lenses unless you reject the dynamically determined es- 
timates of the masses. Including the cosmological uncertainties, the lens data can be fit if 
m* ~ S.llJJ X 10^^ h~^MQ at 95% confidence. This corresponds to a mass-to-light ratio 
of {M/L)b* — 22tl^h at 95% confidence compared to {M/L)b* = (10 ± 2)h in dynamical 
estimates (e.g. van der Marel 1991). The cosmological limits when the masses are high 
enough to fit the observed image separations are nearly identical to the limits we derived in 
the last section using the isothermal sphere model. 

5 The Effects of a Core Radius 

In this section we consider softened isothermal density profiles in more detail to un- 
derstand the effects of a core radius. The three issues we must treat are the dynamical 
normalization of models with a finite core radius, the lensing effects of softened isothermal 
spheres in the presence of magnification bias, and the cosmological effects of introducing a 
finite core radius. 

5.1 The Dynamical Normalization of Softened Isothermal Spheres 

We need to examine the dynamical normalization because the values of (Tdm* and are 
correlated. The core radius reduces the mass near the center of the galaxy and the velocity 
dispersion must increase compared to the value in a singular model to maintain either the 
same observed velocity dispersions or the same average image separations. As a model 
calculation we computed the average line-of-sight velocity dispersion inside one effective 
radius Re assuming a Hernquist (1990) model [p oc r~^(r -|- a)~^) for the distribution of the 
stars with a ~ 0.45i?e- With the assumption that the velocity dispersion tensor is isotropic 
(/3 = in equation (10)), the dark matter dispersion increases as (Tdm oc 1 -|- 2(s/i?e) with 
the addition of a core radius. The numerical coefficient in the scaling law is a function of the 
averaging area, and it rises from 2 to 2.5 if we fix the dispersion inside Re/2 instead of Re. 

Eor a more realistic model we use the 37 galaxies in the van der Marel (1991) sample 
and fit isotropic dynamical models to each galaxy using equations (10) and (11), assuming 
the core radius is a constant fraction of the estimated effective radius for each galaxy, and 
that the velocity dispersion scales as L/ = {(Tdm I c^dm*)'^ ■ The surface of the fit to the 
observed velocity dispersion profiles is shown in Eigure 5, and the dashed line is the scaling 
law estimated from the Hernquist (1990) model. Models with large core radii cannot fit the 
data because of the contradiction between a homogeneous core and a steeply rising luminosity 
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Fig. 5. -(Left) Contours of the for the dynamical fits. The light solid lines show the 68%, 90%, 
95%, and 99% confidence limit changes on Ax^ for the fit to the sample. The dashed line shows 
the expected scaling of (Jdm* with s/Re if we keep the average velocity dispersion interior to Re 
fixed in the Hernquist/softened isothermal sphere dynamical model. 

Fig. 6. -(Right) Variation in cross section (dashed line) and lensing probability (solid line) with 
core radius s for a lens with (Jdm = 250 km s~^. The values are normalized to unity at the 
minimum core radius. The heavy solid line shows the lensing probability excluding image systems 
with detectable central images using the same normalization as for the total probability. The results 
are given for the average over the quasar data sample including selection effects. 

profile. The formal 95% confidence upper limit on the core radius is <^ 0.08i?e or <^ 
300h~^ pc for i?e* = 4:h~^ kpc. While we have greatly oversimplified the dynamical problem 
and ignored the effects of varying the isotropy of the velocity ellipsoid, the basic physical 
effect that the velocity dispersion must rise as the core radius increases is independent of 
any complications. 

For a core radius of = 100h~^ pc the fractional increase in the velocity dispersion is 
7.5% or 17 km s~^, less than the uncertainty in the value of (Tdm*- Nonetheless, its effects 
on models with a finite core radius are striking; it produces a 33% increase in the expected 
number of lenses if we keep the ratio of the core radius to the critical radius fixed {s/b 
constant). This underestimates the effect of increasing (Tdm* because the ratio sjh also 
decreases, leading to a further increase in the number of lenses. As we found in the earlier 
dynamical models of E/SO galaxies in isothermal mass distributions (Kochanek 1992, 1994, 
Breimer & Sanders 1993, Franx 1993), it is absolutely essential either to fit the dynamical 
variables using the observed distribution of image separations or to do the self-consistent 
stellar dynamical problem. 

Self-consistency in the lensing calculation also requires a velocity dispersion that in- 
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creases as the core radius becomes larger. A consistent model of the lenses must keep the 
average image separations fixed as the core radius increases, and the image separation is 
approximately twice the tangential critical radius of the lens (A^ ~ 26(1 — 2/3)). If the core 
radius is increasing (larger /3), then the only way to maintain constant average image separa- 
tions is to also increase the average velocity dispersion (larger b). If we model this by keeping 
the tangential critical radius r_|_ = 6(1 — 2/3) fixed, then the lens cross section r decreases 
as T oc (/3 — 1/2)^ instead of (/3 — 1/2)"^ for /3 ~ 1/2. The condition that the average image 
separation must be constant significantly reduces the decline in the cross section caused by 
the inclusion of a core radius. 

5.2 The Leasing Probabilities for Softened Isothermal Spheres 

Self-consistent calculations of the lensing probability such as Kochanek & Blandford 
(1987), Kochanek (1991), Wallington & Narayan (1992), Kassiola & Kovner (1993), and (in 
most regimes) Maoz & Rix (1993) automatically include the effects of the core radius on 
the magnification bias, but most treatments of softened isothermal models examined only 
the effects of a core radius on lensing cross sections (e.g. Dyer 1984, Blandford & Kochanek 
1987, Hinshaw & Krauss 1987, Krauss & White 1992, Fukugita & Turner 1991, Fukugita 
et al. 1992, Torres & Waga 1995). Core radii have a powerful effect on the cross section 
for multiple imaging and the cross section drops rapidly even when /3 = s/6 is significantly 
smaller than the threshold. Near the threshold of /3 = 1/2, the cross section decreases as 
T oc (1/2 — jSy. However, using the change in the optical depth grossly overestimates the 
effects of a finite core radius on the lensing probability in bright quasar samples. The core 
radius first eliminates images with low total magnifications, but the bright quasar samples are 
dominated by highly magnified images of fainter quasars and magnification bias significantly 
reduces the effects of adding a core radius on the probability. 

We can understand this analytically in the Hinshaw & Krauss (1987) model. The 
tangential critical line of the lens is r_|_ = 6(1 — 2/3)^^^, the radial critical line is r_ = 6(/3 — 
/3V2 - /33/2(4 + /3)i/2)i/2 and the cross section is a^uit = vru^ = TrP{l + 5/3- {3^2 - /S^^^i + 
/3)"^/^/2), where m_ is the radial caustic. Inconsistent models of the effects of a core radius 
estimate the lensing probability by using this optical depth multiplied by the magnification 
bias for the singular model. When images 2 and 3 are merging on the radial caustic, image 
1 is located at ri = Vout = 2/3M_/r^, and the average magnification produced by the lens 
is (Af) = r^^j/u^ = 4/3^/rl. When the core radius is small, the average magnification is 4, 
and as the core radius approaches the threshold of /3 = 1/2 — e the average magnification 
diverges as (Af) oc e~^. The magnification probability distribution is approximately P{> 
M) = [Mo/My for M > Mq where the minimum magnification is Mq = {M)/2 for fold 
caustic statistics (see Schneider, Ehlers, & Falco 1992, Blandford & Narayan 1986, Blandford 
& Kochanek 1987). If we assume a single power law quasar number counts distribution with 
dN/dm oc magnification bias varies with the average magnification as 

B[m) oc (Af)^'^"' for 7 < 0.8. As the core shrinks, the average magnification increases, 
which drives up the magnification bias. The lensing probability, including the change in the 
magnification bias, varies as TB[m) oc e"^"^"' not r oc e"^. For large average magnifications the 
effective value of 7 is the faint slope of the quasar number counts, 7 ~ /3 ~ 0.27 it 0.07 (see 
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Fig. 7.-Cosmological effects of finite core radii. Tfie ieft panei sfiows contours as a function of t fie 
ratio of tfie core radius to tfie effective radius and tfie cosmoiogicai modei. Tfie iigfit soiid fines are 
tfie constraints from fensing aione, tfie dasfied fines adds tfie vefocity dispersion prior as a function 
of core radius estimated in §5.f, and tfie fieavy sofid fines add tfie goodness of fit to tfie dynamicaf 
data. Tfie rigfit panef sfiows dependence of (Jdm* on tfie core radius in an JIq = f cosmofogy. Tfie 
figfit contours use onfy tfie fens data, and tfie dasfied fines incfude tfie prior probabifity distribution 
for (Jdm* estimated from dynamicaf modefs. Tfie contours are drawn at tfie 68%, 90%, 95%, and 
99% confidence fevefs for one parameter. 

§2.5), and TB{m) oc gi-65±0-35_ p^^, ]3j.jg]2t quasars the increase in the bias is greater because 
of the steeper number counts slope, and the effects of the core radius are still smaller. 

To emphasize this point. Figure 6 shows the relative variation of the cross section and 
the true lensing probability including magnification bias for a lens with (Tdm = 250 km s~^ 
as a function of the core radius averaged over the full quasar data sample. For a core radius 
of s = iOO/i~^ pc using the cross section instead of the true probability underestimates 
the lensing probability by about 40%. This comparison still overestimates the effect of a 
finite core radius because we did not include the dependence of (Tdm on s from §5.i. Near 
s = iOO/i~^ pc we estimated that (Tdm should increase by 7%, and the net lensing probability 
with a iOO/i~^ pc core radius may be higher than for a singular model. 

5.3 The Cosmoiogicai Effects of Softened Isothermal Spheres 

We assume that the core radii of galaxies are proportional to their effective radii 
s = s^{L I L^Y'"^ , and the models are characterized by a fixed ratio of s/i?e (as in Fukugita & 
Turner i99i). To simplify the calculations, we set the "Faber- Jackson" exponent to be 7 = 4, 
so the core radius varies with the velocity dispersion as s = s^{(Tdm I c^dm*)'^'^ = ■s^[L/ L^Y''^ . 
We know from §3.i that the effects of the uncertainties in 7 are minimal. If s/s* = 
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{(7dm I (7DM*f^'^ then the strength of the lens varies as /3 = sjh = {s^/b^){aDM / o-dm*T , 
and when u > the core radius suppresses lensing by large [L > L^) galaxies. As we vary 
the velocity dispersion, the maximum tangential critical radius is 

max(r+) = Kfi:^l''v^l\2 + u)-^/^-!/" = 0.4366./3,-°-36 (17) 

when u = 2.8. The maximum image separation is roughly twice the maximum critical radius, 
so to produce images with separations of 2 arcseconds for a source at Zs = 2 and with Oq = 1 
requires cr* ^ 257(s//i~^kpc)°'^"^ km s~^. Thus the scaling law for the core radius introduces 
a sharp threshold on the velocity dispersion as a function of the core radius. 

The right panel of Figure 7 shows the best fit value of (Tdm* as a function of in an 
Oo = 1 cosmology. As expected from §5.1, the velocity dispersion increases as the core radius 
increases with (Tdm* oc 1 + sj R^. This is shallower than the slope seen in the dynamical 
models. The left panel of Figure 7 shows the dependence of the cosmological limits on 
the core radius using only the lens data, the lens data combined with the prior probability 
distribution for (Tdm* derived from the dynamical model of §5.1, and finally the lens data, 
the dynamical velocity dispersion prior, and the likelihood of the dynamical model. If we use 
only the lensing data, the cosmological limits are nearly independent of the core radius - this 
is a radically different picture of the effect of a core radius than that found in inconsistent 
calculations. There is some decrease in the expected number of lenses, and the best fit models 
near Oq = 1 prefer a finite core radius to produce a modest reduction in the expected number 
of lenses. Since the dynamical estimate of the velocity dispersion rises more rapidly than 
the lensing estimate, adding the dynamical estimate of the velocity dispersion eliminates the 
large core radius solutions because of the rapidly increasing number of lenses. If we weight 
the models by the likelihood that they fit the dynamical data, the high of models with 
s/Re ^ 0.05 seen in Figure 5 dominates the limits on the core radius. 

In summary, quantitative estimates of the effects of core radius must include both the 
increases in the velocity dispersion and the increases in the magnification bias caused by 
adding a finite core radius. Most previous calculations have ignored both of these corrections, 
and no previous calculation has included the increase in the velocity dispersion. Figure 7 
covers a much broader range of core radii than is physically reasonable for F/SO galaxies. 
The reasonable range of <^ 200h~^ pc, or s/Re ^ 0.05, is consistent with galaxy dynamics, 
lensed images separations, and the absence of central images (Wallington & Narayan 1993, 
Kassiola & Kovner 1993). Over this range the cosmological limits are independent of the 
core radius. 

6 Spiral Galaxies 

We have hitherto ignored the contribution from spiral galaxies. This is a conservative 
approach, since the additional cross section of the spirals will only strengthen the constraints 
on the cosmological constant. The opposite limit and the maximum contribution is found by 
modeling the spirals as singular isothermal spheres with the characteristic velocity dispersion 
cdm* = (144 lb 10) km s~^ and TuUy-Fisher (1976) exponent 7 = 2.6 it 0.2 used by Fukugita 
& Turner (1991) combined with our more recent estimates of the spiral density = (0.79 it 
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Fig. 8. -Effects of spiraf gafaxies. Tfie sofid (dasfied) fines sfiow tfie fikefifiood contours witfi 
(witfiout) tfie maximaf contribution of spiraf gafaxies as a function of tfie matter density and tfie 
vefocity dispersion of tfie E/SO popufation. Tfie net cfiange in tfie cosmofogicaf fimits is about 
0.05 in fio- Contours are sfiown at tfie 68%, 90%, 95%, and 99% confidence fimits for one variabfe 
in tfie fikefifiood ratio. Tfie fikefifioods incfude tfie fuff uncertainties in tfie number counts and 
"Faber- Jackson" refation for tfie E/SO gafaxies, but onfy tfie number counts uncertainties for tfie 
spiraf gafaxies. 



0.26)h^ X fO~^ Mpc""^ and Scliecliter function exponent a = — f.OOibO.fS from §2.4. Because 
of computationaf fimits we fixed (Tdm* = 144 km s~^ and 7 = 2.6 for the spiraf cafcufation, 
made the variations in the Schechter function exponent for the E/SO and spiraf cfasses fuffy 
correfated (same a for both popufations), and independentfy varied the spiraf density using 
a fog-normaf distribution. There are too many variabfes to independentfy average over the 
prior probabifity distributions for aff parameters of the combined E/SO and spiraf modef 
simuftaneousfy. We varied and Ue independentfy with the broad error bars impiied by 
Marzke et ai. (f994). We couid reduce the variance by using the totai density + Ue and 
the fraction of E/SO gafaxies Up j {n^ + rig) as the independent variabfes since the dominant 
errors come from the fraction of E/SO gafaxies rather than the totai density (see §2.4). At 
the average vaiues of the parameters the modef predicts that 7 ienses are produced by the 
E/SO popufation and f fens is produced by the spiraf popufation when Oq = 1- This is ciose 
to the ratio of the opticai depths given the parameters of the two popufations. Figure 8 
shows that the cosmofogicaf fimits change by onfy AAq — 0.05 when we add the maximum 
contribution from spiraf gafaxies. 

Since even the maximum contribution from the spirais does not significantiy after the 
cosmofogicaf conciusions, we did not pursue their effects in greater detaii for aff cosmofogicaf 
modeis. We did, however, consider how their contribution is modified by extinction and 
finite core radii. We know that spiraf gafaxies contain significant amounts of dust, with 
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central extinctions of order Ab — 1 magnitude (Peletier & Willner 1992) and scale lengths 
of ~ 30 kpc (Zaritsky 1994, Peletier et al. 1995). When we approximate the reddening in 
spirals by a uniform extinction of Ab = 0.5 magnitudes, the expected number of lenses drops 
to one third the number of lenses without extinction. Maoz & Rix (1993) examined spiral 
galaxies with core radii having the same scaling as the effective radii of E/SO galaxies, while 
Krauss & White (1992) examined models with core radii independent of the luminosity. We 
considered the intermediate case where s = s^[cr / cr^y ^ so that /3 = s/b is independent of 
luminosity and all spiral galaxies are equally affected by the finite core radius. We do not 
have to change the values of cr* to compensate for the core radius because the normalization 
is determined by the asymptotic rotation velocity of the galaxies outside the core regions. 
For <^ 70h~^ pc the expected number of lenses changes by less than 20%, but above this 
threshold the numbers begin to drop precipitously. The number of lenses has dropped by 
50% at ~ 200/i-^ pc, and by 90% at ~ bOOh'^ pc. 

Between the effects of the smaller mass compared to E/SO galaxies, reddening, and 
finite core radii, the effects of the spiral galaxies on cosmological limits are negligible. The 
uncertainties in lens statistics must be closer to 20% than 50% before there is any need to 
include spiral galaxies in statistical models. The exception to this rule is for models that use 
large systematic errors to suppress the contribution of the E/SO galaxies by factors of two 
or more; in these models the contribution of the spiral galaxies is important. 

7 Lens Galaxy Luminosity and Extinction 

Lensed quasars differ from unlensed quasars not only because they are multiply imaged 
but because of the luminosity and absorption properties of the lens galaxy. These effects 
are less important for radio selected lenses, but they can introduce large systematic errors 
in statistical analyses of optically selected lenses (Kochanek 1991). In bright quasar surveys 
(m ^ 19 B mags) we do not expect the luminosity of the lens galaxy to be a significant 
problem because the source quasars are almost always significantly brighter than the lens 
galaxy. Reddening, however, can easily reduce the magnification bias by more than a factor 
of two. The question of reddening is whether the reduction in the number of lenses is 
consistent with the properties of E/SO galaxies locally and with other properties of the lens 
systems such as their separations, lens galaxy redshifts, and the similarity of the radio and 
optically selected samples. Eor simplicity we use the nearly singular isothermal models from 
§3. 

7.1 The Luminosity of the Lens Galaxy 

The current optical lens surveys select objects from catalogs of known quasars, so the 
luminosity of the lens galaxy affects the statistical model if it selectively excludes lensed 
quasars from catalogs. At the lens survey stage, the detectability of an extended galaxy 
increases the completeness of the survey because it helps to distinguish the lens from a point 
source. All quasar catalogs have a limit on the luminosity of the lens galaxy beyond which the 
object will not be recognizable as a quasar. Color and variability selected samples are more 
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sensitive to the brightness of the lens galaxy than spectrally selected samples (this includes 
X-ray and radio selected samples). Where identifying the object as a quasar depends only 
on finding quasar emission lines, the galaxy must be significantly brighter than the quasar 
to conceal the strong, broad emission lines. 

The luminosity of the lens is not a serious problem in the current optical lens sur- 
veys because they all selected bright quasars with m ^ 18 B mags. Figure 9 shows 
the approximate effects of the lens luminosity by assigning the lens galaxy the luminos- 
ity Mbj, = —19.9 — 2.57 log ((Tj:)m/cdm* ) and (neglecting k-corrections and evolution for 
simplicity) eliminating all lenses in which the quasar magnitude was less than Am magni- 
tudes brighter than the galaxy. The limit on the brightness of the lens galaxy must exceed 
Am = 3 magnitudes before it introduces any serious systematic bias in the statistical anal- 
ysis, and Am = 3 exceeds any realistic observational limit even for color selected quasars. 
When Ao > the lens galaxy is typically further from the observer and fainter, so for a fixed 
threshold Am you find a larger fraction of the lenses in a high Aq model than in a low Aq 
model. 

We calculated the likelihoods including magnitude limits of Am = 0, 1, 2, and 3 mag- 
nitudes on the luminosity of the lens galaxies. The models included the A;-corrections from 
Coleman, Wu, & Weedman (1980), but no evolutionary corrections. The changes in the cos- 
mological limits are smaller than we can resolve given the coarse gridding of the simulations 
(comparable to the right panel of Figure 3) so there is little point in showing the likelihood 
contours. The Oq = 1 models are the most sensitive to the selection effect; yet when Am = 2 
only a few percent of the lenses are undetectable. Note, however, that the selection effect 
becomes a serious concern in surveys for fainter lenses. For magnitudes between 19 ^ i? ^ 20 
it is an important correction, and for fainter magnitudes [B ^ 20) the statistical results may 
be controlled by the accuracy of the corrections for this systematic bias. This confirms the 
conclusions of Kochanek (1991). 

7.2 Reddening 

Reddening alters the detectability of a lens in three ways. The mean reddening reduces 
the average magnification bias (Kochanek 1991, Tomita 1995), radial gradients in the red- 
dening generally increase fiux ratios (Tomita 1995, Fukugita & Peebles 1995), and patches 
like molecular clouds can absorb single images. In a bright quasar sample the lensed images 
are found near the critical radius of the lens because the magnification and hence the magni- 
fication bias is large there. Since the images lie near the tangential critical radius, gradients 
in the extinction are more important for producing variations in the mean opacity with lens 
redshift than for increasing the fiux ratios. The cross section averaged distance of the critical 
line from the center of the lens galaxy is 2. 0/i~^ ((7/5^/225 km s~^)^ kpc for Oq = 1 and a 
source at Zs = 1. The numerical coefficient rises to 2.7h~^ kpc for Zs = 4. Adding a cos- 
mological constant increases the average impact parameter, with the limiting coefficients of 
3.1h~^ kpc (6.4/i~^ kpc) for Aq = 1 and Zs = I (4). There must be significant dust opacity 
on scales comparable to the effective radius Re to significantly affect lens statistics. Dust 
concentrated in the central regions of the galaxy will have little or no effect. 
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Fig. 9. -The qualitative effects of tlie lens galaxy luminosity (left panel) and extinction (right 
panel) on the expected number offenses as a function of magnitude. The models are normalized to 
estimate the number of lenses per magnitude on the sky normalized by 5 quasars per square degree 
at 19 B mags and a fixed quasar redshift of = 2. In both panels the heavy solid curve at the top 
shows the estimated number counts of unlensed quasars. In the left panel the number counts of 
lenses are shown for JIq = 1 and Aq = 1 with the requirement that the lens galaxy is Am = 0, 1, 2, 
or 3 magnitudes fainter than the quasar. The model includes neither ^-corrections nor evolution. 
In the right panel the number counts of lenses are shown for JIq = 1 and average extinctions of 
Ab = 0, 1 or 2 magnitudes and for Aq = 1 and an average extinction of 2 magnitudes. The model 
assumes the extinction is homogeneous and independent of redshift. The model parameters are 
Ue = 0.0061/i"^ Mpc""^, a = —1, j = 4^ and (Jdm* = 225 km s~^. Only lenses larger than 0.5 arcsec 
in separation and with fiux ratios smaller than 100:1 are included. 

Patchy absorption is the easiest to constrain both in nearby galaxies and in gravitational 
lenses. For example, Goudfrooij et al. (1994) find patchy absorption at the level of Ab ^ 0.2 
mags in their sample of E/SO galaxies. We estimate that the covering fraction / of opaque 
regions in E/SO galaxies must be smaller than / < 0.094 (at 2cr) from the lensed radio point 
sources, because whenever we see one radio image in the optical, we see all the images in the 
optical. We know there is differential reddening in MG0414+0534 (Lawrence et al. 1994) and 
2237+030 (Nadeau et al. 1991), but the absence of detailed spectral and photometric data 
on most lenses makes it impossible to do a more detailed study of the extinction variations 
in known lenses. Direct measurement of the mean absorption is difficult to obtain. In 
spiral galaxies the central reddening in the B band is Ab ~ 1.0 it 0.4 (Peletier & Willner 
1992). The distribution is consistent with Ab oc exp(— r/ro) where Tq ~ 30 kpc both 
by direct observations (Zaritsky 1994, Peletier et al. 1995), indirect estimates from the 
quasar population (Heisler & Ostriker 1988) and damped Lya absorbers (e.g. Lanzetta, 
Wolfe, & Turnshek 1995). The absorption in E/SO galaxies is neither as large nor as well 
studied. Goudfrooij et al. (1994ab, 1995) suggest that there is a diffuse component with 
0.1 ^ Ab ^0.9, but it is unclear whether this an average for the central regions or over many 
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effective radii. Lauer (f988) found no signs of extinction on scales of kpc (limit Ab ^ 0.04 
mag) by modeling the photometry of overlapping cluster ellipticals. The red spectra of 
MG0414+0534 can be produced by 3-6 magnitudes of extinction in the lens galaxy, although 
the uniformity of the extinction for the four images is then remarkable. The red color of some 
radio lenses compared to normal quasars is ambiguous because Webster et al. (1995) find that 
unlensed radio-selected quasars show a broad, fiat color distribution with 2 ^ i? — ^ 10 
(almost as red as MG 0414-|-0534), and O'Dea et al. (1994) find that radio galaxies can also 
show very red colors. Other radio selected lenses show no signs of absorption. In 0957-|-561, 
the radio, B, and R band fiux ratios are all identical even though the A image is n.3h~^ 
kpc from the galaxy and the B image is only 2.6h~^ kpc from the center. This limits the 
differential reddening to AAb ^0.1 mags. 

As a simple model we gave the E/SO galaxies the extinction profile Ab = Ao[l — r/ro) 
for r < To and zero for r > Tq, where Aq is the central value and Tq is a cutoff in the 
dust distribution. Assuming Galactic properties for the dust and gas, a dust to gas mass 
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(eg. Jura et al. 1987, Goudfrooij et al. 1994ab, 1995). Eor a galaxy at distance D with an 
average dust temperature of T = 20 K, the expected 100 fim infrared (Jura 1986, Jura et al. 
1987) and HI radio fiuxes (Knapp, Turner, & Cunniffe 1985; Wardle & Knapp 1986) are 



F[100/im] = O.8A0 
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(19) 

Typical E/SO galaxies have lOO/xm fiuxes (Jura 1986, Jura et al. 1987) and HI fiuxes (Knapp 
et al. 1985, Wardle & Knapp 1986) comparable to these values, implying that the absorption 
in E/SO galaxies locally must satisfy the integral condition that Ao(ro/kpc)^ ^ 1. The 
observation of singular cores in most E/SO galaxies (Tremaine et al. 1994) requires that 
Ao ^ 1, or we would see an apparent core radius at the dust photosphere. Eigure 10 shows 
contours of the dust mass as a function of the central extinction and the scale length, and 
the local evidence suggests that E/SO galaxies occupy the lower left corner. 

To include extinction in the lens calculation we should include the wavelength depen- 
dence of the dust opacity because the redshifted wavelength of an observed B photon was 
further in the UV when it passed through the lens galaxy. EoUowing Heisler & Ostriker 
(1988) and Tomita (1995) we used the Seaton (1979) model for the UV extinction curve 
as a function of wavelength. Converted into a dependence on redshift normalized by the 
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Fig. 10. -Lens probability and average image separation for Aq = 1 and JIq = 1 for tlie extinction 
profile Ab = Ao(l — r/ro) and a 18 B magnitude source at redshift = 2 using the standard E/SO 
lens model. The values of ro and Aq are for an galaxy. The edge radius scales as ro oc (i/i*)^'^ 
and the total dust mass oc AqVq increases as Lj L^. The heavy solid lines show the lines of total 
dust mass equal to lO^M©, lO^M©, 10 ''M© and lO^M©, the light solid lines are the contours where 
the fraction of the lenses found drops to 90%, 80%, 50%, and 10% of the number expected without 
extinction, and the light dashed lines are the contours where the average separation rises to 110%, 
120%, 130%, 140% and 150% of the average separation without extinction. 




Fig. 11. -Likelihood contours for uniform reddening models. The left panel assumes a uniform 
extinction of Ab magnitudes independent of redshift, while the right panel assumes a uniform 
extinction of Ab magnitudes at 2; = and the redshift /wavelength dependence of equation (20). 
The contours show the 68%, 90% confidence, 95%, and 99% confidence limits on one parameter. 
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extinction in the B band at zero redshift, Ab(0), 



Ab{z) = Ab(0) 



1.01 



3.94 + 2.382; + (2.27^-2.33)2+0.280^ 
^.AA -r l.»0^ -r (2.27^-2.33)2+0.280 

10.44 - 4.20^ + 1.53^2] 
- 3.76(/ + z 



z < 0.12 
0.12 < z < 0.61 
0.61 < z < 2.14 
2.14 < z < 3.4 
z > 3.4 



(20) 



where g = \b I 'i ■'iXy For a spectral index of a = —0.5 the effective wavelengths are 
\b = 4403 A and Xy = 5513 A. The fitting formula is comphcated because it models the 
2200 A absorption bump seen in our galaxy, with the peak appearing near z = 1. In this 
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0.5, 1.0, and 1.5 are 1.4, 2.4, and 2.0 times larger than at 



Let nis be the magnitude of the source, then (for two images) the magnitude of the 
lensed images are nii = nis — 2.5 log Mi + Ai and m2 = — 2.5 log M2 + A2 where Mi and 
M2 are the image magnifications with M2 < Mi. We know that the total observed magnitude 
of the quasar is m and the fiux ratio of the images is m2 — nii = 2.5 log M1/M2 + A2 — Ai. 
The two modifications to the lens probability calculation from reddening are to change the 
argument of the quasar number counts ratio in equations (4) and (5) from m + 2.5 log Af 
to nis where is the source magnitude given the magnifications of the images and the 
extinction, 
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and to change the selection effects to use the reddened image fiuxes. 



The simplest model problem is to assume a uniform extinction of Ab at all redshifts and 
see how it affects the magnitude distribution of lensed quasars. Figure 9 shows the qualitative 
effects of extinction on the numbers and magnitudes of lensed quasars. An average reddening 
of 2 magnitudes is needed to reduce the expected number of lenses for Aq = 1 cosmologies to 
the number expected for Oq = 1- The qualitative sign of extinction in the lens population is 
a shift in the break of the lensed quasar number counts to fainter magnitudes. This effect 
cannot be detected in a bright lens survey [B <^ 19 mags) because we do not survey beyond 
the break magnitude, but a fainter lens survey will see significantly more faint lenses in a 
reddened model with a cosmological constant than in a unreddened model with Oq = 1- In 
our example, the difference is a factor of 5 at 20 B mags. 

Variations in the extinction with redshift or radial gradients of the extinction are more 
easily detected than uniform extinction because they modify the distributions of image sep- 
arations or lens galaxy redshifts. For example, models in which E/SO galaxies become nearly 
opaque at higher redshifts (e.g. Fukugita & Peebles 1995) are identical in their lensing ef- 
fects to the galaxy formation models treated by Mao (1991) and Mao & Kochanek (1994), 



®In Heisler & Ostriker (1988) the redshift limit for the first region is incorrectly given as z < 0.19, leading to 
a discontinuous function. 
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allowing them to be ruled out. As a model of the effects of the spatial distribution of dust 
we set the extinction to Ab = Ao[l — r/ro) for r < Tq where the extinction scale length 
varies with the effective radius Tq oc (Z/Z*)^'^, and the total dust mass is proportional to 
the luminosity AqVq oc L. We included the redshift /wavelength dependence of the extinction 
from eqn. (20), and Figure fO shows the effects of dust on the lensing probability and image 
separations for a source with Zs = 2 and m = f8 B mags (the typical survey quasar). A cos- 
mological distribution of lenses with the local dust content of E/SO galaxies {Mdust ^ 10^ Mq 
and Aq < 1) changes the expected number of lenses by less than f 0%. Reducing the number 
of lenses by a factor of two requires M^ust ^ 10"^ Mq, but models at this lower limit alter the 
average image separations by f 0-20% and require very high central opacities. It is impossible 
to reduce the expected number of lenses by a factor of ten without exceeding limits on the 
total amount of metals in the lens galaxy or altering the average image separations by an 
easily detectable amount. 

Because we cannot explore this parameter space fully, we confined our full statistical 
calculation to the effects of uniform extinction. We calculated models both with and without 
the redshift /wavelength dependence of the extinction curve. Remember that these cases are 
the hardest to limit using only the lens data, and they correspond to enormous total dust 
masses - while they lead to the greatest cosmological uncertainty, they are also physically 
unrealistic. We include these two cases only to illustrate the limits of the extinction problem. 
Figure II shows the likelihood contours as a function of the extinction and the cosmological 
model. As expected, there is a degeneracy between the two parameters, since with only five 
bright lenses the magnitude and redshift dependence of the lens probability are too weakly 
measured to constrain the model. The degeneracy is stronger when we neglect the redshift 
dependence of the extinction because it minimizes the side effects of the absorption on the 
image separations. The increase in the opacity with redshift acts like models of galaxy 
evolution where the number of lenses declines with redshift, and these models drive up the 
average lens separation allowing them to be constrained from the lens data (Mao & Kochanek 
1994). If we neglect the dependence of the extinction on redshift, we underestimate the 
amount of extinction required to change the cosmological model by approximately a factor 
of two. 

Both of these models overestimate the extent of the degeneracy. The Aq = I, high 
Ab models are as acceptable as the Oq = I, Ab = models not just because they both 
produce the same number of lenses, but because we put no limits on the reddening in the 
observed lenses. Two magnitudes of extinction in any of the known optical lenses would be 
easily observed, so we know that a significant fraction of the lens population at intermediate 
redshifts is transparent. Moreover, there is no systematic pattern of the radio lenses having 
either higher lens galaxy redshifts or fewer measurable lens redshifts than the optical lenses 
as is required to make the optical and radio selected samples consistent with large amounts 
of extinction and a high cosmological constant. The most powerful constraint on dust in 
F/SO galaxies is probably that quasars with damped Lya absorption systems are not lenses. 
Bartelmann & Loeb (1995) show that there is already a high probability that high column 
density damped Lya systems are lensed by spiral galaxies. If the F/SO galaxies contained 
comparable amounts of gas at 2; ~ I, the probability that the damped Lya systems are 
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lensed would increase dramatically. 



8 The Lens Redshift Test 



The previous calculations did not include information about the relative redshifts of the 
lens galaxy and the source. Kochanek (1992) pointed out that all the known lens redshifts 
were significantly lower than expected for models with a large cosmological constant. The 
analysis neglected both parameter uncertainties and the detectability of the lens galaxy. King 
(1994) and Helbig & Kayser (1995) point out that the selection effects can be important. 

Let pi{z) be the normalized likelihood that the known lens with separation A^^ and 
source redshift Zsi has a lens at redshift z with < z < Zsi^ and let the magnitude of the 
lens at that redshift be mi[z). Let Puii'm-i) be the probability that the lens redshift can 
be measured for a lens of magnitude m in lens system i. Then the probability that we can 
measure the lens redshift for system i is Pm = / dzPMi{'mi)pi[z) and the probability that 
we cannot measure the lens redshift is Pui = 1 — Pm- We can also include the probability 
that the system is found as a lens by adding an additional term to exclude redshifts where 
the lens galaxy is brighter than the quasar, but this is such a negligible part of the total 
probability in the current sample of lenses that we do not include it (see §7.1). 

Suppose we have N lenses with known source redshifts, and that have measured lens 
redshifts and Nu do not. The Kochanek (1992) approach defined the likelihood function by 
Loid = ^f=iPji^j)- Helbig & Kayser (1995) point out that this does not account for where 
the redshifts are measurable, and they advocate only looking at the redshift probability 
distributions over the region in which they are measurable, defining the likelihood function 
to be Lhk = ^^=iPM]{'m)P]{zj) I Pdj- The insensitivity found by Helbig & Kayser (1995) 
is, however, an artifact of the Lhk statistic because it does not account for the fraction of 
the lenses with measurable redshifts. We should be unable to measure the lens redshifts in 
a much higher fraction of the lenses if Aq ~ 1 than if Aq ~ 0. Hence when we include the 
selection effects we should also add the probability of failing to measure the redshifts in the 
systems where we cannot find the lens. We define the likelihood 



which combines the probability of the Njj measured redshifts with the probability that the 
Njj unmeasured redshifts were unmeasurable. By including the information on the fraction 
of the redshifts that are measurable, we restore the sensitivity to the cosmological model. 

We retreat to the SIS model for this analysis. The small core radii permitted by the large 
statistical calculations of §5 will only lead to slight modifications of the tails of the probability 
distribution. The normalized differential probability that a lens has image separation 6 and 
redshift z given a single Schechter (1976) function distribution of galaxies and a "Faber- 
Jackson" relation is (Kochanek 1992) 
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where x = L/ = [AOiDos /'^b^DLsY'^^ is the galaxy luminosity required to produce image 
separation A^^, 6* = 4:Tr[cr£)M*/ ^ and -Dql, Dls and Dqs are the usual proper motion 
distances defined in §2. We average over the uncertainties in a and 7 (as in §3), and we 
include the likelihood of the observed lens separations by adding the term Il^J'-^^^^pk{A.9k) / Pk 
(see §2.3) for the probability that the Nl optical and Nr radio lenses have separations A^^. 
This term is also independent of cosmology (for fiat universes and the SIS model), the galaxy 
number density, and the average reddening. It serves only to constrain the value of (Tdm* 
without needing prior dynamical information. We also compute the likelihood adding the 
standard prior probability for (Tdm*- 

We use the galaxy luminosity model from Kochanek (1992) based on k-corrections from 
the spectral energy distributions of Coleman, Wu, & Weedman (1980) to estimate the R 
magnitude of the lens galaxy as a function of redshift mi[z) (see §2.4). For each lens we assign 
a detection threshold niio (see Table 1). The probability of measuring the lens redshift is unity 
if m < niio and has a Gaussian cutoff with a width of one magnitude exp( — (m — mio)^/2) 
if m > niiQ. The width of the cutoff will compensate for magnitude uncertainties and the 
neglect of galaxy evolution. The correct way to determine the cutoff magnitude is to take 
the best spectrum for each quasar and estimate the threshold at which the lens galaxy 
would be detected, but it is impossible to obtain the raw spectra of many of the lenses. For 
the systems where the lens galaxy is seen (PG 1115+080, 0142-100, MG 1654+1346, and 
CLASS 1608+656) we estimated how much fainter the observed galaxy could be and still 
have a measurable redshift. Depending on the image geometry the limit was either one or 
two magnitudes fainter. For the undetected lenses 1208+1011, H 1413+117, and B 1422+231 
we set the limit to be 2.5 magnitudes fainter than the quasar. For LBQS 1009—0252 and 
CLASS 1600+434 we used stricter limits of R=20 and R=19 magnitudes based on inspection 
of the existing spectra. 

Figure 12 shows the likelihood as a function of Oq and (Tdm*- The Lnew likelihood is 
sensitive to the cosmological model. It sets an upper bound on the cosmological constant of 
Ao < 0.77 at one standard deviation (Aq < 0.90 at 2cr) and the peak is at Aq = 0.4. We expect 
the redshift test to be less sensitive than the full statistical analysis, so these weaker limits 
are not surprising. Large separation lenses have the brightest lens galaxies. In these systems 
the lens redshifts are both known and low (0142-100, PG 1115+080, MG 1654+1346, and 
CLASS 1608+656). Small separation and very bright lenses have lens galaxies too faint to 
measure a redshift in all cosmologies (1208+1011, B 1422+231 and H 1413+117). The two 
intermediate separation systems (LBQS 1009—0252 and CLASS 1600+434) were recently 
discovered. They have not been studied in detail, but the lens galaxies can probably be 
detected in these systems, particularly if Oq ~ 1- As found by King (1994) and Helbig & 
Kayser (1995), the Lhk statistic is very insensitive to the cosmological model, and it cannot 
be used to set any statistically significant limit on Aq. 

9 Conclusions 

The statistics of gravitational lenses, including their numbers redshifts, magnitudes and 
separations are consistent with the expectations for a transparent, constant comoving density 
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Fig. 12. -The right panel shows the likelihood contours of L new for the observed lens redshift sample 
including selection effects. Contours are shown at the 68%, 90%, 94.4%, and 99% confidence limits 
for one variable in the likelihood ratio. The solid lines show the limits using the lensing data alone, 
and the dashed lines show the limits including the prior probability distribution for (Jdm*- The prior 
probability distributions for a and 7 are included in both distributions. The left panel compares 
the cosmological limits for the same selection model using the old Kochanek (1992) statistic, the 
new statistic, and the HK (Helbig & Kayser 1995) statistic. The points on the curves show the one 
standard deviation and 90% confidence limits in the likelihood ratio. Each curve is normalized to 
a peak likelihood of 1. 

of non-evolving E/SO galaxies in a flat Oq = 1 cosmology. The upper limit on the cosmological 
constant in a flat cosmology is Aq ^ 0.65 at 95% confidence. Open universes without a 
cosmological constant and Oq ^ 0.2 have likelihoods near the 90% confidence limit. These 
limits for the first time include all the statistical errors in the input variables determining 
the number and distribution of galaxies and quasars, as well as the statistical uncertainties 
in the lens sample. Is there a cosmic concordance? Perhaps, but not at Aq = 0.65. 

As found by Maoz & Rix (1993) and Kochanek (1995), constant mass-to-light ratio 
models of galaxies modeled with de Vaucouleurs profiles are not consistent with the observed 
separations of the lensed images. An galaxy must have mass 3. or mass- 

to-light ratio [M I L)b* — 22^.1^ h at 95% confidence to fit the lens data even after allowing for 
all the model and cosmological uncertainties. This matches direct estimates of the mass of 
the lens galaxy in MG 1654-|-1346 (Kochanek 1995), but is inconsistent with local dynamical 
estimates (e.g. van der Marel 1991). The cosmological limit found using the de Vaucouleurs 
profile is Aq ^ 0.62 at 95% confidence, little changed from the isothermal models. 

If we model the galaxies as softened isothermal spheres we find that the best fit velocity 
dispersion of the dark matter is (Tdm* — 220 it 20 km s~^ for an galaxy. Unlike the 
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de Vaucouleurs models, we find the same estimates for (Tdm* wlietlier we fit tlie velocity 
dispersions of E/SO galaxies (Koclianek i993, i994, Breimer & Sanders i993, Franx i993) or 
the image separations of gravitational lenses. The best fit models are singular (s* <^ 100h~^ 
pc), consistent with HST photometry of E/SO galaxies (Tremaine et al. 1994), limits on 
the presence of central images in lenses (Wallington & Narayan 1993, Kasiola & Kovner 
1993), models of individual lenses (eg. Kochanek 1995), and general theoretical expectations 
(eg. Dubinski & Carlberg 1991). Adding a core radius has no significant effects on the 
cosmological limits. 

Earlier lens studies were concerned with limits on large cosmological constants (Aq ~ 0.8) 
where the change in the number of lenses is so great that 50% errors in the calculations are 
not very important. Now that we are worrying about Aq ~ 0.5, such inaccuracies are no 
longer tolerable. There are four common problems with lens calculations. (1): The (3/2)^/^ 
correction factor introduced by Turner, Ostriker & Gott (1984) based on simplified models 
of galaxy dynamics does not exist. Models using the (3/2)^/^ correction will overestimate 
the expected number of lenses by 125% and the expected separations by 50%. (2): The 
effects of a core radius cannot be modeled using the changes in the lensing cross section. 
The lensing probability is the product of the cross section and the magnification bias, and 
as the cross section drops the rising magnification bias prevents a dramatic reduction in the 
lensing probability. Eor a core radius of = 100h~^ pc and (Tdm = 250 km s~^, using the 
cross section instead of the true lensing probability overestimates the drop in the probability 
by 40%. (3): In a consistent model, the velocity dispersion must increase when a core radius 
is added to the model. We find the rough scaling that (Tdm oc 1 + ^(s/Re) where 1 < < 3 
depending on the details of the normalization, and s/Re is the ratio of the core radius to 
the effective radius. Because the lens probability is proportional to crj^M^ small increases 
in the velocity dispersion have dramatic effects on the lensing probability. The small rise 
in the velocity dispersion coupled with the effects of magnification bias, cause the weak 
dependence of the cosmological limits on the core radius. (4): The quasar number counts 
model of Eukugita & Turner (1991) uses the slope for the apparent magnitude number counts 
of all z < 2.2 quasars rather than the true luminosity function, causing it to underestimate 
the magnification bias in bright lens surveys. When we use a new number counts model 
derived from Hartwick & Schade (1990), we find that the Eukugita & Turner (1991) model 
underestimates the limits on the cosmological constant by AAq — 0.15. We reach the same 
conclusion with the Boyle et al. (1988, 1990) luminosity function used by Wallington & 
Narayan (1993), or the Boyle et al. (1987) luminosity function used by Maoz & Rix (1993). 

Spiral galaxies are not an important effect in the current lens sample. At most they 
improve the limits on the cosmological constant by about AAq = 0.05. This assumes the 
spiral galaxies have singular cores and no extinction. Given the typical extinction estimated 
by Peletier & Willner (1992) and Zaritsky (1994), their contribution drops by a factor of 
three. Eor core radii larger than ~ 200h~^ pc their contribution drops by a factor of two. 
Between extinction and core radii, there is no need to include spiral galaxies in statistical 
models given the other sources of errors. 

Astronomy is, of course, replete with examples where hidden systematic errors lead to 
incorrect conclusions even though all possible attention was given to the statistical uncer- 
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tainties. We have, however, explored a wide suite of possible systematic errors. Many, such 
as the radial mass distribution of the lens galaxies, ellipticity, core radii, and the luminos- 
ity of the lens galaxy, can be shown to be unimportant in altering the cosmological limits. 
Reddening in E/SO galaxies is a small 10% correction to lens statistics if the extinction is 
limited by local observations of E/SO galaxies to have M^ust ^ 10^ Mq and central opacities 
Ao^ 1. 

The most important remaining class of systematic errors involve evolution, in the num- 
ber, mass, or extinction of the lens galaxies. Some of these issues have already been explored 
(see Mao 1991, Mao & Kochanek 1994, Rix et al. 1994, Eukugita & Peebles 1995), although 
the proliferation of parameters makes a thorough investigation difficult. While it is easy to 
invent a systematic effect that can lower the expected number of lenses and permit a cosmo- 
logical constant, it is hard to do so without altering the observed lens separations, redshifts, 
magnitudes, and the similarities of the radio and optically selected lens samples. A serious 
suggestion for a systematic error must demonstrate consistency with all the lens data. 

As a rule, evolution models that alter the number of lenses by a factor of two will 
inevitably alter some other observable property of the lenses. Eor example, galaxy formation 
or an abrupt rise in the opacity of the E/SO galaxies can reduce the number of lenses only 
at the price of increasing their average separations. This allowed Mao & Kochanek (1994) 
to rule out such evolutionary models as a means of allowing Aq — 0.8. When using evolution 
to avoid the limits on the cosmological constant, we should keep in mind that one effect of a 
large cosmological constant is to make low redshift evolution less plausible (see Carroll, Press, 
& Turner 1992). Moreover, physically plausible merger models tend to alter the distribution 
of image separations but not the expected number of lenses. There is also increasing evidence 
that the E/SO population shows little evolution even at 2; ~ 1 (Lilly et al. 1995). 

The most pressing needs for new observational data (aside from more lenses, and ig- 
noring the difficulty of the observations) are good spectral studies of known lenses, better 
galaxy number counts, and large redshift surveys of the sources in the radio surveys. Good 
spectra or multi-color photometry of the already known lenses, similar to the efforts for 
MG0414-|-0134 (Lawrence et al. 1994), are the simplest way to control the effects of extinc- 
tion. Better spectra are also needed to find the lens and source redshifts, or to set good 
magnitude limits on the presence of a lens galaxy. The number of question marks remaining 
in Table 1 is disheartening. The most important statistical uncertainty in the models is the 
number density of E/SO galaxies, and the lens calculations should really be done using the 
full error correlations between the parameters of the galaxy number counts models (a, n*, 
and L^). Einally, the differential statistics of the radio and optical surveys will be a powerful 
constraint on systematic errors. Yet while the radio surveys are far more productive than 
the optical surveys at finding lenses, the utility of the data is limited by the missing or poor 
characterization of the source population compared to quasars. 
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